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Consistency relations between growth of structure and expansion history observables exist for any 
physical explanation of cosmic acceleration, be it a cosmological constant, scalar field quintessence, 
or a general component of dark energy that is smooth relative to dark matter on small scales. The 
high-quality supernova sample anticipated from an experiment like SNAP and CMB data expected 
from Planck thus make strong predictions for growth and expansion observables that additional 
observations can test and potentially falsify. We perform an MCMC likelihood exploration of the 
strength of these consistency relations based on a complete parametrization of dark energy behavior 
by principal components. For ACDM, future SN and CMB data make percent level predictions for 
growth and expansion observables. For quintessence, many of the predictions are still at a level of 
a few percent with most of the additional freedom coming from curvature and early dark energy. 
While such freedom is limited for quintessence where phantom equations of state are forbidden, 
it is larger in the smooth dark energy class. Nevertheless, even in this general class predictions 
relating growth measurements at different redshifts remain robust, although predictions for the 
instantaneous growth rate do not. Finally, if observations falsify the whole smooth dark energy 
class, new paradigms for cosmic acceleration such as modified gravity or interacting dark matter 
and dark energy would be required. 




I. INTRODUCTION 

A decade after the first firm evidence for the acceler- 
ated expansion of the universe 0, 0] , the study of dark 
energy remains one of the most important yet difficult 
endeavors in theoretical cosmology (e.g. [1, 0, H[). The 
quality of data from a variety of cosmolo gica l probes has 
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to multiple, independent lines of evidence for the accel- 
erating expansion. In the near future, we can expect a 
battery of measurements with unprecedented precision 
that will provide stringent tests of any purported expla- 
nation of cosmic acceleration. 

Despite the tremendous amount of raw information ex- 
pected from upcoming Type la supernova (SN) surveys, 
baryon acoustic oscillations (BAO) from galaxy rcdshift 
surveys, weak lensing, and cluster counting surveys, only 
a handful of parameters associated with the dark energy 
equation of state can be constrained [13, UK . This lim- 
itation arises because most observables depend on cos- 
mological distances and growth, which are integrals over 
the expansion rate, which itself contains an integral over 
the dark energy equation of state. 

Although insensitivity to fine-scale features of the 
equation of state is a drawback for measuring dark en- 
ergy parameters, it is an advantage for testing the con- 
sistency among acceleration observables required by dark 
energy paradigms; since the individual probes of dark 



energy do not depend strongly on the rapidly oscillating 
evolution of the equation of state, neither do the consis- 
tency relations between these observables. For example, 
it is well known that under a cosmological constant ex- 
planation of acceleration or simple paramctrizations of 
the equation of state, distance measurements predict the 
growth of structure in a spatially flat universe. Violation 
of this consistency relation would falsify the standard flat 
ACDM model and its most basic generalizations. 

The goal of this study is to extend these ideas of 
prediction and falsification from simple dark energy 
paramctrizations to general classes of dark energy mod- 
els with time-dependent equations of state, specifically 
scalar field quintessence and dark energy that is spatially 
smooth compared with the dark matter on small scales. 
We start with SN and cosmic microwave background 
(CMB) measurements expected in the next decade and 
make predictions for growth and expansion history mea- 
surements as a function of redshift within dark energy 
model classes parametrized by a complete basis of prin- 
cipal components. Where predictions are tight, observa- 
tions can falsify the model class. Where predictions are 
loose, observations can better pin down the parameters of 
the class, in particular those controlling spatial curvature 
and dark energy that is significant at early times. 

Our study complements previous work [29|, [H HH on 
the observable predictions of classes of dark energy mod- 
els, and follows a long history of studies concerning the 
best way to probe dark energy using cosmological obscr- 
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To our 

knowledge this is the first time that consistency between 
growth, distance, and other expansion observables has 
been studied in such a general and quantitative way. In 
particular, while some previous studies have considered 
predictions of dark ener gy d escribed with a small num- 
ber of parameters (e.g. [29L l3ll . |62|). here we consider 
~ 500 parameters of which about 10 — 15 are necessary 
to completely describe to high accuracy the predictions of 
current or future data. Although we compute equations 
of state for dark energy as an intermediate step between 
distance and growth, we emphasize that our goal is not 
to reconstruct w(z), unlike many previous studies (e.g. 
[Ill [ED, [Zl, H III, M, M, S Izl, Izll)- Rather than 
addressing the quality of constraints on the equation of 
state, here wc are more interested in using dark energy 
parameters only as a tool to study how to falsify basic 
paradigms for cosmic acceleration. 

This paper is organized as follows. In §|TT1 wc describe 
our methods for computing predictions for observables 
from future cosmological data under different dark en- 
ergy paradigms. These predictions can lead to falsifica- 
tion and subsequent generalization of each model class 
as shown in § Mil The main tests of each class are sum- 
marized in § IIVI Appendices provide additional details 
about the inclusion of various data in the construction of 
principal components and in the likelihood analysis, our 
methodology for computing the principal components, 
and tests of the completeness of our paramctrization in 
the growth and expansion observables. 



II. METHODOLOGY 
A. Dark energy principal components 

We parametrize the dark energy equation of state, 
w(z), with a basis of principal components (PCs) [80L [8ll| . 
The PC amplitudes are weighted redshift averages of 
w(z) ordered by how well they are measured, and can be 
straightforwardly computed for a given data set. As we 
discuss in § III Dl and III El principal components based on 
SN distance modulus data in particular provide a nearly 
complete basis for other acceleration observables. We 
therefore treat the PC amplitudes as simply a convenient 
intermediate representation to capture the information in 
the distance modulus and translate it into predictions for 
other observables. 

Principal components can also be defined for other 
redshift-dependent quantities such as the dark energy 
density pue(z) (e.g. [HI, HH) which is related to w(z) 
by 



Pbe(z) = Pci-.o^de exp 



3 / dz 





,1 + w(z') 



(1) 
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where p C r,o is the critical density and Que is the fraction 



of dark energy, both at the present time. Refs. [83j, [84J 
discuss the advantages of using either w(z) or pde(z) 
to describe dark energy. Our choice to use w(z) is mo- 
tivated by the fact that the model classes we consider 
are separated by the allowed values of w: w = — 1 for 
ACDM and — 1 < w < 1 for quintessence. (Values of 
w outside this range are possible for quintessence mod- 
els in which the dark energy density becomes negative, 
but such models are inconsistent with current data as we 
discuss in § 1111 Bl ) 

Specifically, we compute the PCs based on distances 
to Type la supernovae with measurement errors mod- 
eled after the proposed specifications for the Super- 
Nova/ Acceleration Probe (SNAP [85|) experiment. We 
also assume that 300 low-z SNe will be available for cali- 
brating the normalization of the high-z distance-redshift 
relation. Wc supplement the SN observables with con- 
straints on the expansion history at earlier times from the 
CMB acoustic peaks using the precision expected from 
the Planck mission. We take the CMB observables to 
be the matter density scaled to the present, O m /i 2 , and 
the comoving angular diameter distance to last scatter- 
ing, = D{z„) where z* ~ 1090 [24[ . Our assumptions 
about these fiducial SN and CMB data sets and addi- 
tional priors are detailed in Appendix [A] The fiducial 
cosmology we assume for the PC construction (and for 
likelihood analysis) is flat ACDM with present matter 
fraction f2 m = 0.24 and Hubble constant h = 0.73, con- 
sistent with current data [24], HH] . 

The principal component functions &i{zj) are eigenvec- 
tors of the covariance matrix for the equation of state in 
redshift bins {zj}, and they form a basis in which an 
arbitrary function w(zj) may be expressed as 



w(zj) - W Bd (Zj) 



N Z , PC 



(2) 



where on are the PC amplitudes, N z ,pq = 1 + z max /Az 
is the number of bins in redshift, and zj = (j — l)Az. 
We adopt a maximum redshift for variations in w(z) of 
z max = 1.7 to match the largest redshift for our fiducial 
supernova data, and we use a fiducial model Wfid(z) = — 1 
since ACDM is an excellent fit to current data. For these 
choices of z max and Wfid(z), Fig. [T] shows the 15 lowest- 
variance PCs, which form the basis we use for likelihood 
analysis. Wc comment on how the PCs depend on the 
choices of z max and the fiducial cosmology in Appendix iBl 
By normalizing the PCs as 
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the components approach a "continuum limit" as Az — ► 
in which the shapes of all but the worst-determined PCs 
become smooth and independent of Az (or N z pc)- A 
small bin width Az < z^ n also allows us to resolve 
changes in w(z) at rcdshifts below that of the nearest 
supernova in the sample, z^ n = 0.03, which can evade 
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FIG. 1: The first 15 PCs of w(z) (increasing variance from 
bottom to top), with 500 redshift bins between z — and 
Zmax = 1.7. The vertical dashed line shows the minimum 
redshift of the data assumed for computing the PCs, z^ n — 
0.03; SNe are distributed betwee n z^ „ and z msx following a 
SNAP-like distribution (given in [83]) plus a low-z sample at 
z min < 2 < 0.1 (see Appendix |A")| . When computing the PCs 
we include a CMB prior modeled after Planck and marginal- 
ize over fi m , Q,mh 2 , and the absolute magnitude of the SNe. 
The PCs are offset vertically from each other for clarity with 
dotted lines showing the zero point for each component. 



SN constraints as discussed in Appendix [B] We have 
chosen to use bins spaced linearly in redshift, and the 
exact shapes of the components depend somewhat on 
this choice; for example, had we chosen bins with equal 
widths in a = (1 + z)^ 1 or ln(l + z), the weights of the 
PCs would have shifted in redshift [86]. However, the 
most important property of the PCs for our purposes is 
that they form a complete basis for observables such as 
distance and growth (§ III E[) . and this completeness can 
be achieved for a variety of different binning conventions. 

Although we use a large number of redshift bins 
(-Wz,PC ~ 500) to approach the continuum limit of the 
PC shapes, we generally truncate the sum in Eq. © to 
include only the N c < N Zi pq modes that are measured 
best by the fiducial data. Predictions in § IIIII are based 
on a choice of N c = 15, and we explain how this num- 
ber of PCs ensures completeness in various observables 

in §inn 

For model classes that restrict w(z) to some range 
w m ; n < w < w max , we can place priors on the PC am- 
plitudes analogous to those introduced in Ref. [H| for 
reionization principal components. These priors, which 



we define in Appendix [X[ include top-hat bounds on each 
PC amplitude and an upper limit on the sum of squares 
of the amplitudes. Due to the truncation of the num- 
ber of principal components required for likelihood anal- 
ysis, we adopt conservative priors. A combination of PC 
amplitudes is only excluded if the resulting equation of 
state at some redshift exceeds the bounds on w regard- 
less of the amplitudes of the truncated components ({a,} 
with % > N c ). Conversely, satisfying the PC priors does 
not guarantee that a reasonable set of truncated compo- 
nents can bring w(z) back within the bounds. The priors 
therefore include all models within a class, but do not 
necessarily exclude all models outside that class. 

Our baseline dark energy model class is parametrized 
by the PC amplitudes, Q m , and fl m h in a flat universe: 

Abase = {ai, • ■ ■ , oi Na ,£l m , tt m h 2 } . (4) 

The Hubble constant, h = flo/(100 km s _1 Mpc _1 ) = 
(ri m /i 2 /f2 m ) 1 / 2 is a derived parameter in this representa- 
tion. 

In this baseline class we take the dark energy density 
to be constant at z > z max . The underlying assumption 
in this class is that by z = z max the dark energy is al- 
ready much smaller than the matter density as in ACDM, 
and enforcing constant dark energy density at z > z max 
assures that it becomes increasingly irrelevant at higher 
redshift. Note that our baseline model class includes the 
standard ACDM model of a cosmological constant in a 
flat universe, corresponding to {on} = for w$&(z) = — 1. 

B. Early Dark Energy and Curvature 

If observations falsify the baseline model class, we can 
generalize it by including both dark energy that remains 
a substantial fraction of the energy density at z > z max , 
dubbed "early dark energy" [HJlicj, and spatial curva- 
ture. 

To describe early dark energy, we adopt a simple 
paramctrization by assuming a constant equation of 
state, w(z > z max ) — "^oo [861 ] . The dark energy den- 
sity at z > z max can be extrapolated from its value at 

^max aS 

PDE(z) = PDE(Zmax) — " ) • (5) 

\ J- "T ^max / 

This description notably accounts for, but is not lim- 
ited to, scalar field models that "track" at z > z max 
[Oil . l92l , [93^ where the equation of state is determined 
by that of the dominant component, in this case matter 
(w — 0) . We examine the limitations of this paramctriza- 
tion in Appendix [Cl Instead of w^, we use exp(uioo) as 
the parameter for likelihood analysis since models with 
Woo <^ — 1 all have rapidly vanishing dark energy den- 
sity at z > 2 max and are therefore degenerate with each 
other in all observables. We allow the early dark energy 
parameter to vary within the range < exp(w 00 ) < 1, 



3 



where the upper limit eliminates Wrx, > models with 
dark energy density that exceeds the matter density at 
early times. We restrict the allowed range for further 
in model classes where the low redshift equation of state 
is bounded (see Appendix lAl for details). 

Note that components of dark matter that are smooth 
on small scales, for example very light neutrinos, are also 
described by the early dark energy parametrization. We 
will not distinguish between these two possibilities here 
as that would require measurements in a regime where 
cither the neutrinos or the dark energy were not smooth. 

To complete our most general model class, we allow 
for the possibility of spatial curvature, parametrized by 
fix = 1 — £l m ~ ^de- The full parametrization for a dark 
energy model class is therefore 







full 
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,a Nc ,rt m ,fl m h 2 ,cx.p(w 00 ),il K }- (6) 



The present dark energy density Ode is derived from this 
parameter set. Setting Woo = — 1 and Ok = recovers 
the baseline model class of Eq. (|4]) . 



C. Markov Chain Monte Carlo 

We use the Markov Chain Monte Carlo (MCMC) al- 
gorithm to estimate the joint posterior distribution of 
cosmological parameters and derived observables by sam- 
pling the parameter space and evaluating the likelihood 
of each proposed model compared with an assumed data 
set (e.g. see [94| IH, HH). The posterior distribution is 
obtained using Baycs' Theorem, 



V{0\x) 



C{x\0)V{p) 



(7) 



where C{x\0) is the likelihood of the data x given the 
model parameters and V(9) is the prior probability 
density. The MCMC algorithm generates random draws 
from the posterior distribution that are fair samples of 
the likelihood surface. From these samples, we can es- 
timate many properties of the posterior distribution in- 
cluding the mean values, covariance, and confidence in- 
tervals of both the basic set of parameters and derived 
parameters and observables. Convergence of the set of 
random samples to a stationary distribution that approx- 
imates the joint posterior density V{0\x) requires a large 
number of independent samples. We use a minimum of 
four chains per model and determine when these chains 
have a sufficient number of samples for convergence by 
applying a conservative Gelman-Rubin criterion [97j of 
R-l < 0.01. 

The full details of the simulated cosmological data and 
priors used for the MCMC analysis and their likelihood 
functions are given in Appendix [S] and summarized in 

ami 

We assume that all of the fiducial data are consistent 
with a flat ACDM model with fl m = 0.24 and h = 0.73, 
given that this model fits current constraints well. We 



therefore do not consider here the potential for SN and 
CMB data to test this fiducial cosmology. It is, how- 
ever, possible that these future measurements will falsify 
flat ACDM by themselves, even before considering con- 
sistency with additional observables such as growth. We 
have checked that most of our qualitative conclusions do 
not change with allowed alterations of the model under- 
lying the SN and CMB data, and we note exceptions in 
Appendix [Cl 

Given a parametrization for a model class and the fidu- 
cial data, the MCMC posterior distribution then provides 
observable predictions for parameters and derived accel- 
eration observables that can be used as consistency tests 
to attempt to falsify the whole model class. 



D. Acceleration observables 

In this section, we define a set of redshift dependent 
observables that can be probed by future experiments. 
We focus on acceleration observables that can be simply 
computed from the expansion history, leaving for future 
study the detailed relation of these quantities to what 
is actually expected to be measured by specific planned 
experiments. 

We divide the observables into two categories. In the 
first are observables that we assume will be measured by a 
SNAP-like sample of supernovae and the Planck satellite. 
These measurements constitute the data for MCMC like- 
lihood analysis, from which we make predictions for the 
second category of observables in specific model classes. 

Supernova observations constrain the distance modu- 
lus, or relative luminosity distance, between objects of 
different redshift in the sample. We take the SN data as 
a starting point since of the known methods for constrain- 
ing the acceleration, it has the finest resolution in redshift 
and hence its principal components form the most com- 
plete set for providing testable predictions. Supernovae 
have an additional advantage of being sensitive to low 
redshifts z < 0.5 where dark energy dominates the en- 
ergy budget and where other probes like BAO and weak 
lensing do not have enough volume and distance, respec- 
tively, in order to strongly constrain dark energy. We 
take the same Planck CMB constraints on £l m h 2 and an- 
gular diameter distance as used in the PC construc- 
tion. 

The SN and CMB data make predictions within a 
model class for the remaining observables, which include 
the expansion rate H(z), the absolute distance D{z), the 
growth function G(z), and the growth rate f(z). 

The expansion rate, allowing for a general dark energy 
component and spatial curvature, is 



H(z) = H Q 



3 , PDE(z) 



Pcr,0 



1/2 



(8) 

Except when dealing with CMB observables, we generally 
ignore the contribution of radiation to the expansion rate 
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since it is a negligible fraction of the density at low z. 
The absolute distance observable we use is the comoving 
(angular diameter) distance 



D(z) 



(|ttK|H 2 )V2 



Sk 



2x1/2 



dz' 



, (9) 



where the function Sk (x) is equal to a; in a flat universe 
(J7k = 0), sinhec in an open universe (f2x > 0), and sin a; 
in a closed universe (fix < 0). Luminosity distances, 
whose ratios are measured by SNe (see Appendix [X)) . are 
simply related to Eq. © by d h (z) = (1 + z)D(z). 

We define the growth function as G(z) oc (1 + z)Di(z), 
where D\(z) = 8(z) / 8(z- m it) describes growth of the 
matter overdensity 8 normalized at an initial rcdshift 
-Zinit = 1000 during matter domination. If £l m (z) = 
Q m (l + z) 3 {Hq/H(z)} 2 = 1, then D x (z) oc (1 + z)- 1 and 
G(z) is constant. The growth function obeys 



G" 



G' 



;tt m (z) 



G = 0, (10) 



where primes denote derivatives with respect to In a. We 
normalize G to its value at z- m i t , taking the initial condi- 
tions to be 



--(1 - Woo)^DE(zinit), (11) 
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G(zinit) — 1, G'(zinit) 



where G' follows from a power-law solution to Eq. (fl"0|) 
assuming riDE(^init) <C 1 and neglecting curvature since 
it has little effect on the expansion rate at early times. 
We do not include radiation when solving Eq. (fit))) , and 
we assume that the dark energy component is smooth on 
scales below the horizon at z < Zmit- 

Because growth measurements can be compared either 
to low redshift data sets to obtain the relative growth 
Gq{z) = G(z)/G(z = 0) between a redshift z and the 
present, or to recombination through the CMB acoustic 
peaks to obtain G(z), we show predictions for both in the 
following sections. The latter will ultimately be limited 
by the measurement of the optical depth to rcionization 
from CMB polarization due to the translation between 
the observed acoustic peak amplitude and the intrinsic 
fluctuations at recombination. 

The logarithmic growth rate is defined as 



_ dlngi G' 



(12) 



which is commonly approximated as f(z) = QZ, (z) where 
the growth index is 7 « 0.55 for flat ACDM [ill, [§g[. 
Measurements of 7 have been proposed as a way to test 
general relativity; we examine this idea in the context of 
various classes of cosmological models in § MI Dl 

In this paper, we remain agnostic about the techniques 
that best probe these observables and simply assess the 
precision to which they can be predicted in certain model 
classes. However, some caveats are useful to keep in 
mind. Although we make predictions as a function of rcd- 
shift that can include fine-scale features, measurements 
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FIG. 2: Redshift dependent quantities for the fiducial flat 
ACDM cosmology with f2 m = 0.24 and h — 0.73 as assumed 
for PC construction and for the default data sets for MCMC, 
including the fractions of the total density in matter and dark 
energy, Q m (z) and $1de(z) (top; solid blue and dashed red, re- 
spectively), comoving angular diameter distance D(z) (mid- 
dle, solid blue), inverse of the expansion rate H~ 1 (z) (middle, 
dashed red), growth function relative to early times (bottom, 
solid blue), and growth rate / = 1 + dlnG/dlna (bottom, 
dashed red). The vertical dotted line is plotted at 2 max = 1.7, 
the maximum redshift for the PCs and for SNe in the likeli- 
hood analysis. 



will typically only constrain coarse-grained averages of 
the predictions over wide bands in rcdshift. For exam- 
ple, observations of the imprint of baryon acoustic oscilla- 
tions (BAO) on galaxy clustering in directions transverse 
to the line of sight provide a standard ruler to constrain 
absolute distances D(z), and BAO measurements along 
the line of sight can constrain the expansion rate, H(z). 
However, a volume of > 1 Gpc 3 is required to obtain 
accurate measurements of cither quantity, resulting in 
smearing in redshift. 

Likewise, weak lensing measures both the growth func- 
tion and ratios of distances, but the broadness of the 
lensing kernel and scatter in photometric redshifts again 
prevents a purely local measurement. Growth rate mea- 
surements that involve the redshift space distortion of 
galaxy surveys and galaxy bias information from lensing 
suffer from broadening from both data sets. Growth and 
growth rate measurements from the cluster abundance 
or weak lensing in the nonlinear regime also involve inte- 
grals over the past history o f gro wth and not merely the 
instantaneous linear growth [lOOj. 

Finally, it is useful to place weak current priors on 
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parameters related to the observables. In the broadest 
model classes that we consider, the SN and CMB mea- 
surements alone are not sufficiently predictive to elimi- 
nate even highly deviant models. To ensure that we only 
include models that arc not already ruled out by observa- 
tions, we include priors on the fractio n of d ark energy at 
recombination f2DE(.z*) from WMAP [ 1 1 j ] . the absolute 
distance D{z = 0.35) from current BAO measurements 
from SDSS [l3[, a nd th e Hubble constant Ho from the 
HST Key Project |l02| . The former represents the im- 
pact of the change in the growth function near recombi- 
nation on the first few acoustic peaks of the CMB. We 
conservatively do not include the expected improvement 
on this measurement from Planck. The BAO distance 
prior serves mainly to reduce the possible deviations from 
a flat geometry for dynamical dark energy models with 
spatial curvature. The role of the Hubble constant prior 
in this analysis is to limit the variation in w(z) at very 
low redshifts as we describe in Appendix [B] 

In summary, the observables that we predict are the 
expansion rate H(z), comoving absolute distances D(z), 
the growth history G(z) relative to recombination or 
G = G(z)/G(z = 0) relative to the present, and the 
growth rate f(z). The redshift evolution of these quanti- 
ties for the fiducial flat ACDM model is plotted in Fig. [2] 
Predictions from SN and CMB data are derived by con- 
structing these observables from models in the MCMC 
samples described in § III C[ and we explore the implica- 
tions for various dark energy model classes in § IIIII 



E. Completeness 

To make predictions that can reliably be used to falsify 
paradigms for dark energy, our parametrization must de- 
scribe any effects that models within the class might have 
on observables. In particular, we must ensure that the 
set of principal components that parametrize variation in 
w(z) form a complete basis for representing changes in 
growth and expansion observables relative to the fiducial 
cosmology due to changes in the dark energy equation 
of state. In this section we summarize our criteria for 
completeness, and refer the reader to more detailed dis- 
cussions of these issues in Appendices [B] and ICl 

When computing the PCs from the Fisher matrix for 
SN and CMB distances, we have a choice of whether to 
fix or marginalize over the parameters other than the 
binned w(z). These decisions affect PC shapes due to de- 
generacies between w(z) and the other parameters. We 
choose to fix fix and w^, thereby including degenera- 
cies with curvature and early dark energy among the 
well-measured PCs. We marginalize J7 m in the PC con- 
struction, thus reducing the completeness for represent- 
ing sharp transitions in w(z) at z < z m ^ n as we describe 
in Appendix [Bl Since these transitions are largely indis- 
tinguishable and limited mainly by the external Hubble 
constant prior, completeness is not important here. We 
therefore choose in this instance to sacrifice completeness 



for efficiency in representing the well constrained redshift 
range. 

Given the set of PCs, the next question is how many 
out of the full set of N Z) pc components we need to keep 
as parameters (see § 111 A[) . Note that neglecting even the 
high variance PCs can have large effects on the equation 
of state. However, since all of the observables contain 
integrals over w(z) the effects of these rapidly oscillat- 
ing PCs (sec Fig. [TJ tend to cancel out for the redshift 
dependent quantities of interest, especially the distance 
and growth observables. In general, the number of com- 
ponents necessary for completeness, N c , will be larger 
than the number of dark energy parameters that can be 
measured to some specified accuracy. We typically find 
that the predicted range of observables changes fraction- 
ally by less than a few percent between MCMC analyses 
with 10 and 15 PCs (see Appendix EJ). The agreement 
is somewhat worse for H(z) in some cases, but discrep- 
ancies occur mostly at z < z max where oscillations of 
H(z) about the fiducial model would be averaged out in 
BAO measurements over wide bins in redshift. We con- 
clude that N c « 10 is sufficient for completeness, but we 
present results from the larger set of 15 PCs to further 
reduce any remaining artifacts due to incompleteness. 

The completeness of the parametrization could in prin- 
ciple depend on the choice of fiducial model that we adopt 
to represent the true cosmology for future observations. 
We assume that the future SN and CMB data will be 
consistent with wad = —1 as is true of current measure- 
ments, and we examine an alternate choice of fiducial 
model in Appendix [Cl 

Completeness may also depend on the redshift range 
over which the PCs are defined, < z < z max . The 
choice of z max influences the definition of early dark en- 
ergy as well, since we ascribe any deviations from w = — 1 
behavior at z > z max to early dark energy parametrized 
by constant w = Woo. We find that choosing z max = 1.7 
to match the redshift coverage of the fiducial SN data 
nicely balances between defining the PCs where cosmo- 
logical data have significant support (which argues for 
a lower z max ) and having the PCs be a complete rep- 
resentation for other observables (arguing for a higher 
Zmax)- Additionally, our parameter is not intended 
to be a complete description of early dark energy but 
rather a means of monitoring its observable signatures. 
We present tests of both our choice of z max and the early 
dark energy parametrization in Appendix [Ul 



III. TESTING DARK ENERGY PARADIGMS 

Distance measurements from SNe and constraints from 
the CMB make predictions for the acceleration observ- 
ables described in § III Dl that can be tested by future 
experiments. These predictions are made within the con- 
text of a paradigm for acceleration, e.g. a cosmological 
constant. Where the predictions are weak, the observ- 
ables can be used to estimate parameters within the class, 
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FIG. 3: Illustration of the model classes (large circles) and 
subclasses (represented by points within different regions of 
the circles) that we consider in this paper. Open points mark 
the initial (simplest) type of model within each class, and 
arrows indicate paths to the more complex models in the 
class. The ACDM class does not contain models with signifi- 
cant early dark energy since the dark energy fraction vanishes 
rapidly at high z for all allowed cosmological constant models. 



and where they are strong, precision measurements can 
potentially falsify the whole dark energy paradigm. 

In the following sections, we step through the predic- 
tions for various dark energy model classes. Guided by 
Occam's razor and criteria for falsifiability, we begin with 
the simplest model that satisfies current constraints: flat 
ACDM. Models of this type make the firmest predictions 
and are therefore easiest to falsify. The next simplest 
and most predictive model class is ACDM with one ad- 
ditional parameter, spatial curvature. This more general 
class is particularly interesting in that falsification would 
rule out a cosmological constant. 

Since ACDM corresponds to a constant dark energy 
equation of state w = — 1, our next step in generalizing 
the class of models is to allow w to vary with redshift. 
We first consider a restricted range, — 1 < w < 1, corre- 
sponding to the allowed values of the equation of state in 
quintessence models where a canonical scalar fiel d is r e- 
sponsible for cosmic acceleration at late times (e.g. jl03l |). 
Within the quintessence class, we study how the predic- 
tions change when early dark energy and nonzero curva- 
ture are added to the basic model. 

Instead of going directly from a cosmological constant 
to arbitrary equations of state, one could test several 
intermediate models along the way such as constant uu 



or the tw o-parameter model w(a) = wq + w a (l — a) 
jl04 Il05| . and many authors have used this approach 
in analyzing cosmological data. Here, however, we take 
the view that unless the dark energy falls into a well de- 
fined physical classification such as a cosmological con- 
stant or a canonical scalar field, there is no reason to 
define a particular functional form of w(z) as a class and 
not another. 

For the final model class, we allow w to vary over a 
much wider range than in the case of quintessence but 
retain the requirement that dark energy is smooth com- 
pared with dark matter on scales associated with the 
measurements of growth. For example, non-canonical ki- 
netic terms can lead to equations of state with w < — 1 
jl06j | . If such a field has a sound speed substantially be- 
low the speed of light, then the growth predictions pre- 
sented here would only apply below its sound horizon 
fT07j . 

We allow an arbitrary but large range of the equa- 
tion of state within Aw — 4 of the fiducial w = — 1, 
so —5 < w < 3. We limit the range of w to enable the 
MCMC sampling to converge to the joint posterior distri- 
bution of the parameters more easily. This range is large 
enough to include extreme departures from ACDM and 
quintessence models, particularly considering the conser- 
vative nature of our priors on PC amplitudes (see Ap- 
pendix [X]). As with the quintessence model class, for the 
more general class of smooth dark energy models we also 
examine the effects of early dark energy and curvature 
on predictions for obscrvablcs. 

Figure [3] shows a Venn diagram representation of the 
model classes that we study. The ACDM model class 
forms the innermost circle since it is a subset of all the 
other dark energy classes. Quintessence occupies a larger 
portion of the model space since the equation of state is 
allowed to vary within the range — 1 < w < 1, and models 
with even more general equations of state are contained 
within the outer "smooth dark energy" circle. The four 
quadrants of this diagram separate models that are ei- 
ther flat or have nonzero curvature and that either do or 
do not have a significant fraction of early dark energy. 
Within each of the three model classes, an open circle 
marks the simplest type of model (flat with no early dark 
energy), and arrows lead to other models within the class 
with additional degrees of freedom. Since the require- 
ment that w = —1 at all times for ACDM models implies 
that dark energy is always negligible at high redshift (see 
Fig. [2]) , the ACDM class has only two types of models in- 
stead of four. See § IIVI for an index relating these model 
classes to the figures in the following sections. 

The plots we present in the following sections show 
the range of fractional deviations in observables relative 
to the fiducial flat ACDM model that is allowed by the 
assumed SN and CMB data, based on the distribution of 
MCMC samples. We plot these predictions for various 
growth and expansion observables (§ III D[) at rcdshifts 
< z < 4. Any future measurements that fall outside 
the predicted range of values would falsify a model class. 
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Since it is impractical to show the whole posterior dis- 
tribution at a number of redshifts for several different 
observables, we plot only the 68% and 95% limits of the 
distributions and a single example model selected from 
the MCMC samples. The confidence limits are defined so 
that the probability (i.e. number of samples) is equal at 
the upper and lower limits, with 68% (or 95%) of the sam- 
ples between those limits. This definition corres pond s to 
the "minimum credible interval" (MCI) of Ref. [l08j]. A 
useful property of these MCI limits is that the confidence 
region includes the mode of the samples even when the 
distribution is strongly skewed. 

The redshift range in our plots of observable predic- 
tions extends beyond the coverage of the assumed SN 
sample so that we can make predictions for observables 
at higher redshift. Such high redshift observations are 
especially important for limiting the effects of curvature 
and early dark energy parameters that can change ob- 
servables at redshifts beyond the reach of the SN data 
set, although the CMB distance prior can also play this 
role in simpler classes of models. 

There are two ways in which a model class may still be 
falsified even if future growth and expansion observations 
appear to be consistent with the SN and CMB predic- 
tions for that class. Observations that point to redshift 
evolution of an observable that is inconsistent with the 
evolution in the majority of samples would exclude the 
corresponding model class despite appearing to be con- 
sistent with the predictions at any single redshift. Simi- 
larly, measurements of multiple observables could falsify 
a model class if they are inconsistent with the predicted 
correlations between those observables. Since these types 
of inconsistencies between predictions and observations 
can be difficult to see in the types of plots shown here, 
we use the example models in each figure to help point 
out some of the trends in redshift and between observ- 
ables. 

In this study we do not address the feasibility of mak- 
ing the measurements required to falsify various model 
classes using future data sets. Instead, we focus here 
on determining what types of observables are the most 
effective at distinguishing competing theories for accel- 
eration and what kind of precision in their measurement 
would be required. We leave the task of connecting our 
results with realistic expectations for upcoming dark en- 
ergy probes for future work. 



A. Testing ACDM 

The flat ACDM model has only two free parameters, 
f2 m and Ho, whose values are closely tied together by 
the CMB prior on fl m h 2 . With this simple model, the 
fiducial SN and CMB data, or even the CMB data alone, 
make strong predictions for the other observables (see 
Fig. [p. 

The uncertainty in the growth function G(z) at red- 
shifts approaching recombination is zero by definition, 



and only increases to 0.5% by z — (quoting 68% CL 
here and throughout this section). The expansion ob- 
servables at z max , D(z = 1.7) and H(z = 1.7), are pre- 
dicted with 0.4% and 0.2% accuracy, respectively. At 
low z, D and H have equal fractional uncertainties (since 
lim 2 ^o D(z) = z/ Hp), corresponding to an accuracy of 
0.7% for Hq [l09j | . In comparison, current estimates of 
Hq have uncertainties of 3.8% from WMAP alone and 
1.9% from combined WMAP, SN, and BAO measure- 
ments [24J. Note that there is an extremely tight and 
potentially falsifiable prediction for H at z ~ 1 of 0.09% 
in flat ACDM. This prediction is driven mainly by the 
tight CMB distance prior which effectively reduces the 
remaining freedom in ACDM to one parameter. In a 
flat universe, D = J dz/H(z), so allowed variations in 
H 2 oc fl m h 2 (l + z) 3 at high z must be compensated at 
low z by an opposing variation in Hq to preserve . 

Given these strong predictions, which are in large part 
already available from current data, any future detec- 
tion of deviations in growth, absolute distance, or the 
expansion rate at the percent level at any redshift would 
provide evidence against flat ACDM. The predictions 
are driven mainly by the CMB prior and hence the SN 
data themselves can be viewed as a stringent test of 
flat ACDM (see e.g. [ill). Conversely, testing the flat 
ACDM predictions on other observables does not depend 
strongly on having the SNAP SN data set in hand. 

Future observations that rule out flat ACDM would in- 
dicate a need for additional complexity in the model. Al- 
though the predictions in Fig.Ushow what measurements 
would falsify flat ACDM, they do not indicate what kind 
of generalizations of the model would give such alternate 
predictions. For example, tight predictions of AD/D at 
z > 3 may not be so interesting if there are no reasonable 
models that can generate deviations from the flat ACDM 
predictions there. Even in the context of falsifying flat 
ACDM it is important then to look at the predictions of 
an extended class of models. Of the possible directions 
for generalizing the flat ACDM model outlined in Fig. [31 
we first examine the effects of including spatial curva- 
ture. We then consider the alternate option of allowing 
time variation of the equation of state in the following 
sections. 

For ACDM with curvature, growth and expansion pre- 
dictions from the fiducial data are roughly a factor of 2 
weaker than for flat ACDM but remain at the percent 
level (see Fig. [4}. The impact of the SN data is much 
greater when we allow curvature to vary, since the CMB 
constraints alone are no longer sufficient to fix the ob- 
servables at low redshifts. The maximum uncertainty, at 
z = 0, is 1.1% for G and 1.7% for H from D and H. 
The pivot point in H at z ~ 1 also disappears, leading to 
the largest fractional change in the precision of predic- 
tions. In terms of measuring or constraining curvature 
under the ACDM paradigm, the redshifts with the weak- 
est predictions and the largest change from flat ACDM 
offer the most fruitful epochs for measurement. Con- 
versely, the redshifts with the strongest predictions offer 



8 



o 

o 
<1 



o 

\ 

o 

o 
< 



Q 



\ 

■x 




0.02 



_ 1 1 1 1 1 1 1 1 
- n K *o 


i | i i i i | i i i i _ 




















o 

o 
<1 



o 

\ 

o 

o 



Q 
< 



\ 
<] 




0.04 

: -l^wgi 

0.02 r ACDM 




FIG. 4: Forecasted predictions from SNAP SN and Planck 
CMB data for growth and expansion observables, showing the 
influence of curvature on predictions for ACDM models. The 
growth function G is defined relative to its value at recombina- 
tion, and Go is defined relative to the present value. Absolute 
distance D differs from SN relative distances due to uncer- 
tainties in Hq. Shaded regions enclose 68% CL regions and 
curves without shading are upper and lower 95% CL limits, 
plotted as fractional differences from the fiducial flat ACDM 
cosmology. The model classes are flat ACDM (light gray) and 
ACDM with nonzero spatial curvature (dark blue). An exam- 
ple model with nonzero curvature is also shown (dashed red 
curve). Figures [5T- 1 1 1 1 and 1 141 all follow the same format. 



the best opportunity to falsify the cosmological constant 
altogether. As with flat ACDM, with current con strai nts 
these predictions weaken only by a factor of ~ 2 . 

In summary, under the assumption of ACDM, with or 
without curvature, SN and CMB observations make firm 
predictions <~ 1% for growth and expansion observables 
at all redshifts. Future measurements that rule out a cos- 
mological constant as the source of cosmic acceleration 
would represent a significant advance from a standpoint 
of fundamental physics. To proceed beyond this point 
and determine the best observables and redshifts to tar- 
get in order to distinguish among alternate models for 
acceleration, we need to widen the model class again. 



FIG. 5: Effects of generalization of flat ACDM (light gray) 
to quintessence (dark blue; example model: dashed red). 
Quintessence is defined to have — I < w < 1 with w(z) 
parametrized by 15 PCs. Here Woo — —1 and Q.k — to 
eliminate early dark energy and curvature. 



B. Testing quintessence 

If measurements of growth or expansion observables 
exclude ACDM as a viable model, the remaining dark en- 
ergy model classes are ones where the dark energy equa- 
tion of state at late times (z < z max ) is a free function of 
redshift (see Fig. [3]). We add this freedom to the models 
by parametrizing w(z < z max ) with PCs as described in 
§ IH A[ We use the first 15 PCs of w(z) at z < z max in 
the MCMC likelihood analysis as this number suffices for 
completeness in the observables to percent level precision 
(see 5 III El and Appendix EJ). 

As long as the scalar field potential remains positive, 
the equation of state for quintessence is bounded in the 
interval — 1 < w < 1. Negative potentials that violate 
this bound in the past either would not produce the re- 
quired acceleration or would display easily falsifiable fea- 
tures. Our implementation of the canonical scalar field 
prior as described in Appendix [X] is very conservative; 
all quintessence models are allowed by the prior but not 
all models allowed by the prior can be represented as a 
canonical scalar field. 
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We begin in Fig.[5]with the predictions for quintessence 
in a flat universe with Woo = — 1 to eliminate early dark 
energy. Ruling out flatness in the context of ACDM 
does not necessarily mean that curvature is required in a 
more general model class, so we begin with the simplest 
quintessence models from Fig. Sand generalize to models 
with curvature and early dark energy later. 

There are several notable features of quintessence pre- 
dictions when compared with those in the ACDM class 
of models. The limits on growth are no longer monotonic 
with rcdshift and in particular show that growth suppres- 
sion at the 1 — 2% level at z ~ 1, which would rule out 
ACDM, is allowed in the more general quintessence con- 
text. Distances are typically predicted at a level that is 
about twice that of flat ACDM at z < z max and compara- 
ble to ACDM with curvature. The exception is at z = 
where the quintessence class allows for sharp changes in 
the equation of state at z < z^ n as described in Ap- 
pendix [Bj The prior on the Hubble constant restricts the 
amplitude of such changes. Conversely, with a prior on 
the quintessence model class excluding these sharp tran- 
sitions, precision Hubble constant measurements would 
play the same role as absolute distance D(z) or H(z) 
measurements at z^ n < z < 0.1 here and in all of the 
following cases. 

Finally, the quintessence class allows a substantially 
wider range of Hubble parameter predictions H(z). How- 
ever, most of the allowed variation comes from rapid os- 
cillations in the dark energy density, demonstrated by 
the example model in Fig. [5l As mentioned in § III Dt 
observations such as BAO that constrain H are spread 
over a wide rcdshift bin, averaging out much of this oscil- 
latory behavior. Constant shifts in the average AH/H ~ 
—AD/D that would be allowed by the SN measurements 
are still highly constrained in this model context by the 
CMB distance measurement. Without curvature or early 
dark energy, the absolute distance between z max and re- 
combination is nearly fixed and thus the Planck measure- 
ment also fixes shifts in the distance scale below z max (see 
Appendix [B| . 

The flat quintessence class of models with no early dark 
energy on the whole remains highly predictive and falsifi- 
able. We next examine what kind of observations might 
falsify this class by requiring early dark energy or non-flat 
geometries. 

Fig. [H] shows the predictions for quintessence in a flat 
universe with Woo ^ — 1, i.e. with early dark energy. Re- 
markably, the predictions for absolute distance obscrv- 
ables remain nearly unchanged. In particular, there is 
no substantial increase at z > z m ax nor are compensat- 
ing H and D shifts at z < z max allowed. The lack of 
additional freedom in the observables is mainly due to 
restricting the equation of state to — 1 < w < 1. Early 
dark energy allows increased H at z > z max and hence 
decreased distance between z max and recombination. To 
remain consistent with the CMB measurement of -D*, 
these changes must be compensated by a reduction in 
the average AH/H and an increase in absolute distances 
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FIG. 6: Effects of early dark energy [dark blue; exam- 
ple model: dashed red) on quintessence models [light gray). 
Quintessence models from Fig. [5] are generalized to have Woo 
vary, with S1k = to eliminate spatial curvature. 



below z max . Given that the fiducial model is ACDM with 
w = — 1, a roughly constant negative shift in H requires 
dark energy to decrease with redshift, i.e. a "phantom" 
equation of state with w < — 1. Since such values of w 
are not allowed in the quintessence class, predictions for 
D and H remain tight even with early dark energy. 

On the other hand, although growth predictions are 
still at the ~ 1 — 2% level they show an interesting fea- 
ture that is a signature of early dark energy. Growth 



suppression at z 



> 



is allowed at a larger level and 



results in a nearly constant offset in growth at lower red- 
shifts (see dashed curve in Fig. [6]). Positive AG/G is not 
allowed since the amount of early dark energy can only 
increase from the fiducial model, which has almost no 
early dark energy due to the assumption of Woo = ~ 1- 
Growth relative to z = 0, AGq/Gq, is largely unaffected 
by the extra freedom allowed by early dark energy. The 
smoking gun of early dark energy is therefore a compo- 
nent to the growth function deviation that is nearly flat 
at z ^ z max . 

Fig. [7] shows the predictions for quintessence in a non- 
flat universe with Woo = — 1, i.e. with no early dark en- 
ergy. Here the first difference is the change in the rel- 
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FIG. 7: Effects of curvature (dark blue; example model: 
dashed red) on quintessence models (light gray). Quintessence 
models from Fig. [5] are generalized to have £Ik vary, with 
Woo — — 1 to eliminate early dark energy. 



FIG. 8: Effects of curvature and early dark energy (dark blue; 
example model: dashed red) on quintessence models (light 
gray). Quintessence models from Fig. [5] are generalized to 
have both f?K and Woo vary. 



ative growth Go between z max and z = that is now 
allowed to be several percent. While both early dark en- 
ergy and curvature can suppress growth relative to the 
early matter dominated epoch, this shift in Go is a unique 
signature of nonzero curvature. A measurement that in- 
dicates AGo/Go ~ 2 — 4% at z ~ z max would falsify fiat 
quintessence models, with or without early dark energy. 
A measurement beyond this level (or with the opposite 
sign) would falsify non-flat cases as well. 

The second main difference due to curvature is that 
high rcdshift negative deviations in D(z) at z > z max are 
now allowed at the ~ 4% level. A small curvature affects 
the distance to recombination more than it does distances 
at lower redshift. In an open universe, becomes larger 
and therefore allows dark energy the freedom to compen- 
sate by introducing a constant shift down in AD/D and 
up in AH/H at z < z max . The allowed amplitude of the 
shift is limited by our BAO prior at z = 0.35. Note that 
since the shift is constant in redshift, any single absolute 
distance measurement at z < z max also suffices to con- 
strain this mode. An example of such a nearly constant 
shift is shown by the dashed curve in Fig. [Jj 

In a closed universe, the required compensation at z < 



Zmax to preserve -D* is in the opposite direction, and as 
with early dark energy the w > — 1 quintessence prior 
limits this possibility. Thus one-sided deviations in D, 
average H, and growth relative to z = are signatures 
of curvature in the quintessence model class. 

Fig.[5]shows the predictions for quintessence with both 
curvature and early dark energy. Here a much greater 
range of early dark energy densities is allowed since cur- 
vature in an open universe can compensate for the de- 
creased distance to recombination due to early dark en- 
ergy. The main difference is a large increase in the al- 
lowed nearly constant offset in the growth G relative 
to recombination at z < z max which can now approach 
20%. The growth deviations are in fact limited by our 
early dark energy prior from WMAP (see Appendix [A"| : 
without this prior the growth offset could have reached 
~ 40%. Large suppression of the growth relative to re- 
combination indicates both early dark energy and curva- 
ture in the quintessence context. 

Note that the growth at z < z max relative to z = re- 
mains nearly as well predicted as in Fig. [7J with no early 
dark energy and hence is equally falsifiable. Likewise, 
predictions for D do not weaken further because constant 
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deviations in distance are limited by the BAO prior on 
D(z = 0.35) in both cases. On the other hand, allowed 
models with significant early dark energy do weaken pre- 
dictions for H at z > z max . 

Even in the most general quintessence class, there 
are still a few firm, percent level predictions. Neither 
the growth G nor the distance D can be appreciably 
larger than the ACDM prediction, although both can 
be smaller. Lower average H than in the fiducial flat 
ACDM model is not allowed at z < z max due to the 
w > — 1 bound and at z > z max due to the CMB prior 
on Q, m h 2 . Suppression of growth relative to high z at a 
level of > 5% at z = z max must remain nearly constant 
at z < z max . Observations that violate these predictions 
would falsify the entire quintessence model class. 



C. Testing smooth dark energy 

Falsification of quintessence would challenge many the- 
ories of dark energy and motivate consideration of more 
complicated models than single canonical scalar fields. 
Our generalization to the smooth dark energy class en- 
compasses equations of state with — 5 < w < 3 and 
requires that the dark energy remain smooth compared 
with the matter on scales associated with growth mea- 
surements. 

Figures l9T-fTl"l and [14] show the growth and expansion 
predictions from SN and CMB data relative to the fidu- 
cial model for the class of general smooth dark energy 
models. As with quintessence, we present predictions for 
models both with and without curvature and/or early 
dark energy. Remarkably, Figure [9] shows that for flat 
models without early dark energy the effect of dropping 
the quintessence bounds on w(z) weakens predictions by 
less than a factor of two. Thus the more general class 
of smooth dark energy without curvature or early dark 
energy is nearly as falsifiable as flat quintessence. 

Including early dark energy in smooth dark energy 
models relaxes the predictions in ways qualitatively sim- 
ilar to the quintessence case, but with allowed devia- 
tions that are somewhat larger at 68% CL and noticeably 
larger at 95% CL. The additional freedom in growth and 
distances comes mainly from dropping the lower bound 
on w(z) and allowing phantom dark energy models with 
w < —1. This change enables a reduction of the dark 
energy density at z < z max from the fiducial w = — 1 
model, which can compensate for and thereby allow a 
higher fraction of early dark energy while maintaining 
consistency with the CMB distance prior. Since dark en- 
ergy is the dominant component at late times, decreasing 
its density typically results in lower total density at low z, 
which increases absolute distances there. The extra early 
dark energy suppresses the growth at z > z max , and the 
reduced dark energy density at low z raises growth back 
up toward the fiducial model slightly by z = 0. Com- 
paring Figs. [5] and [TUl we find that these models with 
significantly lower H(z < z max ) are the only flat smooth 
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FIG. 9: Effects of generalizing — 1 < w < 1 quintessence 
models (light gray) to smooth dark energy with —5 < w < 3 
(dark blue; example model: dashed red). The smooth dark 
energy model class generalizes the quintessence models from 
Fig. [5] using the same 15 PCs for w(z < z max ) with Woo = 
— 1 and fix = to eliminate early dark energy and spatial 
curvature. 



dark energy models that would not already be excluded 
if observations had falsified all quintessence models. 

In contrast, allowing spatial curvature to vary in the 
smooth dark energy model class produces models with 
qualitatively new types of behavior not seen in the 
quintessence models. These new observable deviations 
from flat ACDM provide smoking-gun signatures of cur- 
vature and an equation of state beyond quintessence. 

Although some of the model classes examined so far 
allow growth relative to early times to be suppressed by 
~ 20% or more compared with the fiducial flat ACDM 
model, none of the previous cases allow growth to be 
enhanced by more than ~ 2%. However, models with 
weak bounds on w(z) and nonzero spatial curvature can 
have G(z) be as much as 10 — 15% higher than the fiducial 
model at z = 0, as shown in Figure 1111 Likewise, the 
growth Go relative to z = can be lower than in the 
fiducial model by 5 — 10% or more at z > z max , whereas 
it was previously limited to deviations of at most a few 
percent in this direction. Both effects correspond to an 
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FIG. 10: Effects of early dark energy (dark blue; exam- 
ple model: dashed red) on smooth dark energy models (light 
gray). Smooth dark energy models from Fig.[9]are generalized 
to have w^, vary, with £7k = to eliminate spatial curvature. 



FIG. 11: Effects of curvature (dark blue; example model: 
dashed red) on smooth dark energy models (light gray). 
Smooth dark energy models from Fig. [9] are generalized to 
have Qk vary, with io~o = — 1 to eliminate early dark energy. 



enhancement of growth at z < z m ax that occurs in closed 
models (Ok < 0). 

The reason why such models become viable when we 
abandon the — 1 < w < 1 prior is similar to the expla- 
nation for the differences in the early dark energy pre- 
dictions of quintessence and smooth dark energy mod- 
els. Closed universes have smaller distances to recom- 
bination, so without some other means to increase the 
distance, the models with Ok < are inconsistent with 
the CMB distance prior for the fiducial cosmology Re- 
moving the quintessence bounds on w(z) allows for lower 
dark energy density (with w < —1) at low z, which in- 
creases the total distance to last scattering and allows 
closed models to match the CMB constraints. The lower 
dark energy density at z < z max and resulting enhance- 
ment of distances are also reflected in the predictions for 
H and D in Fig.[TlJ Having Ok < and Ode {% < z max _) 
lower than the fiducial value means that O m (z <C 2 max J is 
higher than in the fiducial model, which also contributes 
to the additional growth of structure at low redshift. 

Although the average H(z) at low redshift can be re- 
duced considerably relative to flat ACDM, predictions 
for H(z) at high redshift have a sharp lower limit. The 



CMB constraint on O m /i 2 places a strong lower bound on 
H(z) corresponding to the expansion rate in an Einstein- 
de Sitter universe. Closed models can exceed this bound 
slightly since the curvature slows the expansion, but for 
allowed values of Ok this is a small effect. This lower 
limit is present in all previous model classes as well but 
its impact is less visible. 

Like the non-flat quintessence predictions of Fig. [JJ the 
predictions for smooth dark energy models with curva- 
ture are asymmetric about the fiducial model. In fact, as 
Fig. [11] shows, the confidence regions are so skewed to- 
ward the closed models that the fiducial model lies on the 
edge of the 68% regions. This is in spite of the fact that 
the fiducial model has the maximum likelihood by defini- 
tion, and that the definition of confidence limits we use 
is chosen to include the peak probability. Moreover, the 
allowed regions of the observables in Fig. [TTJ are skewed 
in the opposite direction of the predictions for the cor- 
responding quintessence models in Fig. [JJ and therefore 
appear inconsistent with those predictions, despite the 
fact that the quintessence class is a subset of the more 
general smooth dark energy class. The reason for these 
discrepancies is that predictions in the non-flat cases with 
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FIG. 12: Effects of priors on the degeneracy between curva- 
ture and w(z) PCs for the smooth dark energy models with 
curvature (but not early dark energy) in Fig. 1111 Top: ID 
marginalized posterior probability P(Qk) for the default pri- 
ors that are flat in PC amplitudes on (solid, dark blue), and 
for alternate priors that are flat in the density of each PC at 
Zmax relative to z = 0, pi(z m ax)/pi(0) as defined in Eq. (fl3)l 
(dot-dashed, red). The dotted curve is the mean likelihood 
distribution for flat Qi priors. Bottom: Probability contours 
of 51k vs. ct-2 at 68% and 95% CL for the same priors on PCs 
as in the top panel. Dashed lines mark the fiducial values, 
fix = and ct2 — 0. 
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FIG. 13: Effects of priors on smooth dark energy models with 
nonzero curvature. Dark blue: flat top-hat priors on oa as in 
Fig. 1111 Light gray: flat priors on the density of individual 

PCs, /5i(z max )/pi(0). 



large variations in w(z) both above and below w = —1 
are so weak that the shapes of priors on the PC ampli- 
tudes become important in determining the extent of the 
confidence regions. 

The influence of priors on the predictions is illustrated 
in Fig. 1121 where we show distributions of S1k for two dif- 
ferent choices of priors. With our usual top-hat prior on 
{cti}, the posterior probability for Ok is strongly skewed 
toward closed models. The distribution of the mean like- 
lihood of MCMC samples, on the other hand, is peaked 
at the fiducial value of £!k = as expected (dotted curve 
in Fig. US]). 

The discrepancy between posterior probability and 
mean likelihood can be tr aced to a volume effect in the 
parameter space (e.g. see |ll2j |). Models with more neg- 
ative J7k have a wider range of values of cc* that fit the 
data well; this is demonstrated by the banana-shaped 
contours for S7k and «2 in the lower panel of Fig.[T2J and 
other PCs show similar widening of the parameter vol- 
ume at more negative J7k- When other parameters are 
marginalized to obtain the ID posterior distribution for 
f^K, or for one of the growth and expansion observables 



at some redshift, the result is a skewed distribution. 

The basic reason for this volume effect is that the dark 
energy density depends exponentially on u>(z), which is 
a linear combination of the PCs, so changes in {at} at 
small pde have less effect on observables than changes at 
large pde- To test how much the observable predictions 
are affected by the priors, we use alternate priors that 
are flat in the contribution of each principal component 
to the dark energy density at z max relative to z = (see 
Appendix IX)) . 



Pi (^max) 

pW) 



exp 



3a, 



dz 



di(z) 
l + z 



(13) 



so that the total dark energy density at z max is 
) = Pde(O) Yl l [Pi{z max )/p l (Q)}. Figure US shows 
that with this new prior, the predictions become more 
symmetric around the fiducial model and also allow mod- 
els that are acceptable under the quintessence subclass. 

These ~ la shifts indicate that the predictions from 
cosmological data alone are so weak that the exact con- 
fidence region depends on arbitrary theoretical priors on 
the measure in dark energy model space. Therefore, any 
measurement of these observables at a level of precision 
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FIG. 14: Effects of curvature and early dark energy (dark 
blue; example model: dashed red) on smooth dark energy 
models (light gray). Smooth dark energy models from Fig. [9] 
are generalized to have both Woo and Qk vary. 



FIG. 15: Effects of priors on smooth dark energy models 
with curvature and early dark energy. Same as Fig. 1131 but 
comparing priors for the dark blue model in Fig. 1141 



comparable to the predictions is interesting in the con- 
text of smooth dark energy with curvature regardless of 
the sign of the deviation from flat ACDM. At the same 
time, conclusive falsification of general smooth dark en- 
ergy models with curvature would require much larger 
deviations where there are no models with good likeli- 
hood values. 

Note that the dependence on PC priors is only sig- 
nificant for classes of models that allow w < — 1 and 
have nonzero curvature. For all of the previous cases - 
ACDM, quintessence, and fix = smooth dark energy - 
the confidence limits of observables shift by only < 1% 
when we switch from one set of PC priors to the other 
(except for H, for which the limits change by up to a 
few percent at some redshifts near z max ). The volume 
effect is not a consequence of our particular parametriza- 
tion of w(z); for example, there is a similar shift toward 
JIk < for non-flat dark energy models parametrized as 
w(a ) = wp + w a (l — a) when the priors are flat in wq and 

w a [mi. 

Using the priors that arc flat in the PC amplitudes, the 
addition of early dark energy to smooth dark energy mod- 
els with curvature (Fig. [1~4"|) appears to make little differ- 



ence to the qualitative predictions, except in some of the 
95% limits. However, much of the effect of early dark 
energy on these models is masked by the flat-{ai} PC 
prior. As Fig. [TSI shows, changing the priors to be flat in 
Pi(zmax) / Pi(0) affects the observable predictions in ways 
that are similar to the previous case without early dark 
energy, e.g. making many of the distributions of allowed 
models more symmetric around the fiducial model. With 
this alternate prior we see that the additional freedom in 
early dark energy combined with nonzero curvature en- 
ables models with significant growth suppression relative 
to high z and increased H(z > z max ) to fit the assumed 
data sets, as it does in the corresponding quintessence 
predictions of Fig. 

The sensitivity to the priors and the weakness of the 
predictions in general means that despite the great poten- 
tial of future observations for measuring spatial curvature 
and early dark energy in models with general equation of 
state variation at low redshift, statements about falsifi- 
cation of the entire smooth dark energy model class must 
be made with care. 
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D. Beyond Smooth Dark Energy 

One lesson from the analysis in the previous section is 
that the combination of general variation in the dark en- 
ergy equation of state with early dark energy and nonzero 
spatial curvature allow a wide variety of cosmological 
models to fit future SN and CMB data. Falsification of 
the most general smooth dark energy model class there- 
fore appears to be quite difficult, especially given the 
dependence of the growth and expansion observable pre- 
dictions on priors that must be set arbitrarily in the ab- 
sence of an underlying theory for dark energy. However, 
even these very general models make some firm predic- 
tions about the relations between observables that could 
potentially be falsified by future measurements. 

One example of a robust prediction for the observ- 
ables is the following: given the flat ACDM model that 
matches the SN and CMB data and assuming that the 
dark energy always contributes positively to energy den- 
sity, H(z) at z > 2 can be no more than ~ 5% lower than 
in ACDM. This limit was noted in § 1111 C|) in the context 
of non-flat smooth dark energy models with = — 1, 
and it still holds when we include early dark energy. 

More interestingly, the predicted redshift evolution of 
expansion and growth observables still exhibits certain 
regularities. Notice in Fig. [JJJ that before we introduce 
curvature, the growth history, absolute distances, and 
expansion rate are predicted at the few percent level by 
SN and CMB data. Freedom in the spatial curvature 
greatly reduces the precision of these predictions. How- 
ever, the curvature is set by a single parameter J7k with 
well defined effects on each of the observables. By taking 
advantage of our knowledge of the impact of curvature, 
we can effectively regain much of the predictive power 
that exists for flat models. 

As an example, consider the growth function. In the 
most general model class, the majority of the freedom 
in growth comes from curvature and early dark energy. 
We can distinguish between the two by noting that the 
deviations in growth from the fiducial flat ACDM model 
have different redshift dependence, as illustrated by the 
sample growth histories allowed by either curvature or 
early dark energy plotted in Fig. [16J The main effect of 
early dark energy is to suppress growth by a constant fac- 
tor at early times, so growth functions in a flat universe 
with varying amounts of early dark energy have similar 
shapes but different amplitudes at low z. On the other 
hand, curvature tends to have a more gradual effect on 
growth continuing to z = 0, and also allows more en- 
hanced growth relative to the fiducial model than early 
dark energy. 

An observed growth history that cannot be described 
by some combination of the effects of curvature and early 
dark energy would present a major challenge to the dark 
energy paradigm. For example, growth at z < z max rela- 
tive to today (Go) that is > 5% higher than expected in 
fiat ACDM would be difficult to explain with dark energy 
for allowed values of the spatial curvature. 
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FIG. 16: Growth functions of MCMC samples in the smooth 
dark energy model class (— 5 < w < 3) that include either 
early dark energy (EDE) at z > z ma x (w x 7^ —1; top), cur- 
vature (SIk 7^ 0; middle), or both (bottom). We plot growth 
relative to early times in the left column of panels, and growth 
relative to the present on the right. Dashed red curves show 
growth in the fiducial flat ACDM model. Samples are selected 
randomly from those with likelihoods satisfying A% 2 < 4, 
but for visual clarity we plot samples that are approximately 
evenly spaced in G(z — 0) (left) or Go(z — 4) (right). The 
dotted vertical line in each panel marks 2 = z max . 



Moreover, deviations in Gq(z) are nearly a one param- 
eter family that is ordered by curvature f2j<- The top 
panel of Fig. [T7] shows the correlation between J7k and a 
linear combination of Go(z max ) and G(z max ) in smooth 
dark energy models. The curvature Qk mainly depends 
on Go(zmax), but the degeneracy in D* between AH/H 
at z < z max and early dark energy introduces a small 
correction since such changes in the expansion rate affect 
the growth rate at low z and therefore change Go- The 
G(z max ) term can correct for this degeneracy since the 
amplitude of G is sensitive to the amount of early dark 
energy. 

By comparing the combination of Go and G correlated 
with curvature at multiple redshifts, one can test the gen- 
eral class of smooth dark energy models with early dark 
energy and curvature. For example, the lower panel of 
Fig. [T7] shows these combinations of growth observables 
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FIG. 17: Comparison of S7k and linear combinations of the 
growth function relative to flat ACDM, AGo(z) and AG(z), 
at 2 = 2 max and 2 = 1 for randomly selected models (with 
A\ 2 < 4) in the smooth dark energy class with early dark 
energy and curvature. The coefficients used here are K\(z = 
1) = 0.24, K 2 (z = 1) = 0.10, Ki{z = 2max ) = 0.17, and 
K 2 (z — 2 max ) = 0.09. Note that these specific values may 
not produce accurate estimates of SIk for fiducial cosmologies 
other than the one assumed here. 



at z max = 1.7 and at z = 1; observations that give dif- 
ferent values for the linear combinations of Go and G at 
these two redshifts would falsify this most general model 
class. 

This regularity in the growth relative to today can 
be viewed as a generalization of tests involving the 
linear growth rate f(z). In particular, the relation- 
ship f(z) = f22a(z) with 7 0.55 has been proposed 
as a potentia l test of a ll smooth dark energy models 
@ Ell Gil EE E3- In Figs. HHH2ni we plot the 
predictions from the forecasted SNAP supernovae and 
Planck CMB data for the growth rate and growth in- 
dex for selected model classes from the previous sections, 
allowing a redshift dependent growth index 



7 (z) = 



MM*)]' 



(14) 




Note that the growth index is somewhat different from 
our other obscrvables since measurement of 7 requires not 
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FIG. 18: Effects of curvature on the ACDM growth rate. 
Predictions from future SN and CMB data for the growth 
rate / = 1 + dlnG/dlna, plotted relative to the fiducial 
model, and for the growth index 7 = In // ln[f2 m (z)]. The 
model classes here are ACDM either assuming flat geometry 
(light gray) or allowing nonzero curvature (dark blue; exam- 
ple model: dashed red). Shaded regions enclose 68% of the 
allowed models and curves without shading are upper and 
lower 95% confidence limits (not visible for flat ACDM pre- 
dictions of 7 due to the tightness of those constraints). The 
three example models plotted here and in Figs. 1 191 and 1201 arc 
the same as the ones in Figs. [4] [8] and 1141 respectively. 



only the growth rate but also some method for determin- 
ing the fractional matter density Q m (z) at the same red- 
shift. A measurement of H(z) combined with the CMB 
constraint on fl m h 2 could provide an estimate of the lat- 
ter quantity, since fl m (z) oc O m /i 2 (l + z) 3 H~ 2 (z). 

In the context of ACDM, Figure [TS] shows that both 
/ and 7 are precisely predicted by future SN and CMB 
data as expected given the tight constraints on other ob- 
scrvables within this model class. The growth index is 
nearly constant, with a small slope at low redshifts when 
the cosmological constant dominates, and deviations of 
~ 0.5% by z = 4 for models with nonzero curvature. 

In the context of quintessence, both predictions weaken 
substantially as shown in Fig. 1191 The growth rate is 
not as well predicted as the difference in growth between 
3 = and z max used in Fig. [T7] Like H , it effectively has 
only one integral over the time-varying equation of state 
instead of two as for distances and the integrated growth 
history. It is therefore equally sensitive to features in 
w(z). 

In fact, the redshift dependence of the predictions for 
f(z) in Fig. Q33 closely mimic those for H(z) in Fig. [SI but 
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FIG. 19: Effects of curvature and early dark energy on 
the quintessence growth rate. Same as Fig. [18] but for fiat 
quintessence models with Woa = — 1, i.e. no early dark energy 
(light gray), and quintessence with both nonzero curvature 
and early dark energy (dark blue; example model: dashed 
red). 



with opposite sign. Since £l m (z) is tied to H(z) via the 
CMB prior on Ct m h , the similarity between deviations 
in / and in H suggests a strong connection between / 
and Q m (z), consistent with a constant value of 7. How- 
ever, the predicted values of 7 shown in the lower panel 
of Fig. ll9l cover a much wider range for quintessence than 
for ACDM, with significant variation of 7 with redshift 
in some models. For flat models without early dark en- 
ergy, the extra freedom in 7 only appears at z < z max , 
but including curvature and early dark energy allows 5% 
deviations in 7 above the fiducial value at z > z max , only 
slightly less than the uncertainty in / or H at these red- 
shifts. These nearly constant deviations in 7 at high z are 
well approximated by the ex pect ed dependence on early 
dark energy given by Ref. [114j . modified for nonzero 
curvature, 

7(z > z max ) rj — — / DE + -/k, (15) 

5 — bwoo 7 

where /j = Qi(z)/[l — Q m (z)]. For 17k = and Woo = 
— 1 we recover the usual growth index from this formula, 
7 = 6/11 w 0.55. 

One of the more interesting features of the predictions 
for 7 in the quintessence class is the widening in the un- 
certainty at z <~ 1. As Eq. (fTS")) shows, deviation from 
w = — 1 changes the value of 7 at high redshift. Likewise, 
low-redshift variation in w can perturb the growth index 



FIG. 20: Effects of curvature and early dark energy on the 
smooth dark energy growth rate. Same as Fig. [19] but for 
smooth dark energy models (—5 < w < 3). 



from its usual value of 7 w 0.55. At z < 0.5, SN distance 
constraints ensure that while w may oscillate rapidly with 
redshift, it never deviates far from w — — 1 on average. 
However, as z approaches z max the constraints from SNe 
weaken, allowing w to vary significantly from —1 over 
longer periods of time; this variation enables 7 to deviate 
from the ACDM prediction. Furthermore, at z > 1 where 
fl m (z) is typically near unity, ln[r2 m (z)] is close to zero 
and therefore the value of 7 derived from Eq. (|14[) is more 
sensitive to small changes in the relation between f(z) 
and il m (z). For example, the model plotted in Figs. [5] 
and 023 has significant changes in H and / between z = 1 
and z = 1.5 relative to the fiducial flat ACDM cosmology, 
leading to a bump in 7 for this model at z ~ 1. 

Generalizing to smooth dark energy models with — 5 < 
w < 3 reveals another potential problem with using 7 to 
test smooth dark energy. Recall from the previous section 
that a closed universe is allowed for models in this class, 
and is favored if priors are flat in the PC amplitudes. In 
closed models where J7k is sufficiently negative, fl m (z) 
can cross unity with Q m (z) > 1 at high z and Q m (z) < 1 
at low z, and the same is true of f(z). Since these two 
functions do not cross at exactly the same redshift due 
to the slight lag between density and growth, the growth 
index of Eq. p4[) has a singularity when £l m (z) = 1. 
The dashed curve in Fig. [201 is one example of such a 
model. As a result, the predictions for 7, particularly 
in the tails of the distribution, blow up at z > 1 where 
these singularities occur. This effect is an artifact of the 
7 parametrization since 7 can take any value when / w 1 



18 



and fl m (z) w 1, but it makes it difficult to interpret limits 
on 7 beyond z ~ 1 for the most general class of smooth 
dark energy models. 

Although constant 7 remains a good approximation for 
many dark energy models, large variations in the dark en- 
ergy density or spatial curvature with f2j< < can weaken 
the link between f(z) and fl m (z) (or H(z)). An observed 
deviation from the expected value and near constancy of 
7 would certainly falsify ACDM and some simple dark 
energy models, but using the growth index as a test of 
smooth dark energy in general may require refinement 
of the standard parametrization of Eq. (fT4|) to account 
for significant w(z) variation and crossing of Q m (z) = 1. 
Alternatively, one can adopt the more general approach 
of examining the integrated growth function at various 
redshifts as discussed above. 

There are also other, more fundamental but more qual- 
itative means of testing smooth dark energy. By defi- 
nition, on scales where the dark energy remains smooth 
there is no particular scale for growth in the linear regime. 
Models of acceleration that involve coupling of dark en- 
ergy to dark matter or modifications of gravity that in- 
troduce new scales in addition to the Hubble scale gene r- 
ically imply scale dependent linear growth (e.g. [1171 ]). 
Such models can also feature differences in dynamical 
and lensing mass measurements. Gcncrically, modified 
gravity models that satisfy local constraints on gravity 
also break the relationship between the linear and non- 
linear growth of structure again by the i ntroduction of a 
new scale to the problem jll8l . I119L Il20l | . 

Another way in which the standard cosmological 
paradigm might be falsified is through observed viola- 
tions of the relation d^(z) = (1 + z)D(z) between lu- 
minosity distances and comoving angular diameter dis- 
tances. Examples of mechanisms for violating this re- 
lation include photon-axion mixing, photon decay , and 
non zero torsion in the gravity theory (e.g. jl2ll . Il22l 
Il23| ). Therefore, the "duality relation" between the two 
distances is an interesting test of exotic new physics pos- 
sibly related to acceleration. There also exists a more 
general (but related) consistency relation between the 
comoving distance D(z) and the Hubble parameter H(z) 
that holds in any homogeneous and isotropic FRW model 
jl24| and can be tested using accurate cosmological ob- 
servations of the two functions at any redshift. Our stan- 
dard assumptions could also be falsified through observed 
violations of homogeneity or isotropy signaling a break- 
down in the validity of the FRW metric, or by observing 
time variation in fundamental constants. 

We have shown that dark energy degrees of freedom 
and spatial curvature permit the basic distance, expan- 
sion, and growth observables to vary greatly. Neverthe- 
less, there are still many ways in which the dark energy 
paradigm for acceleration could be falsified. 



IV. CONCLUSIONS 

Using a combination of quantitative tools including 
principal components of the dark energy equation of 
state and MCMC analysis with simulated future data, 
we demonstrate that combined constraints on dark en- 
ergy from the measurement of distances, growth, and 
the expansion rate provide many ways to test not only 
specific dark energy models, but also general classes of 
models. In particular, a high-quality supernova sample 
such as that anticipated from SNAP and the CMB data 
expected from Planck make strong predictions for other 
observables in the context of a wide variety of models. 
Follow-up observations of the predicted observables of- 
fer the exciting possibility of falsifying model classes and 
demonstrating the need for a new paradigm for acceler- 
ation. 

Figure [5l] summarizes the basic dark energy model 
classes and their generalizations, arranged by the allowed 
values of w and the total number of MCMC parameters 
[Eq. ©]. In general, observations that falsify classes to- 
ward the upper left corner of this "tree" of model classes 
require adding more freedom to the models by moving 
down and/or to the right in the tree. This diagram serves 
as an index to the figures in § Mil each line between a 
pair of model classes is labeled by the number of the fig- 
ure in which a comparison of observable predictions for 
those two classes may be found. 

In the context of the current standard model, flat 
ACDM, predictions for acceleration observables come 
mainly from the CMB. The SNAP SN data themselves 
provide one stringent test of flat ACDM predictions from 
Planck. Other types of observations are also well suited 
for testing flat ACDM. In particular, constraints from fu- 
ture SN and CMB data have a narrow pivot point in the 
expansion rate where H(z = 1) is predicted to ~ 0.1% 
(at 68% CL), making such a measurement an especially 
interesting target for future BAO experiments and other 
probes of H. Likewise, growth observables are predicted 
to better than 0.5% at all redshifts. 

Even if we drop the assumption of flatness, ACDM 
remains highly predictive. The pivot point in H(z) pre- 
dictions disappears when nonzero spatial curvature is al- 
lowed, but the ACDM model class can still be falsified 
with 1% deviations in any of the acceleration observ- 
ables at any redshift. 

Remarkably, allowing general time variation of the 
dark energy equation of state w(z) within the class of 
quintessence models (— 1 < w < 1) does not significantly 
weaken predictions if we assume flatness. For these mod- 
els, SN and CMB data predict growth and absolute dis- 
tance to ~ 1 — 2% precision. 

Including the possibility of early dark energy, e.g. a 
scalar field that tracks the matter density at high red- 
shift, does not weaken these predictions substantially as 
long as the SN and CMB observations remain consistent 
with a cosmological constant as current data suggest. An 
increased fraction of early dark energy reduces the dis- 
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Number of MCMC parameters 




FIG. 21: Index for dark energy model classes compared in figures in § Mil The number of parameters varied in the MCMC 
likelihood analysis increases from left to right. For extensions to the baseline model within each class, the additional parameters 
(Ok, Woo, or both) are listed. Red numbers along lines connecting two models indicate the number of the figure in this paper 
in which we plot growth and expansion predictions for that pair of models. 



tance to recombination, Z?», and must be compensated 
by either allowing nonzero curvature or reducing the dark 
energy density at late times. The latter option requires 
w < — 1 which is forbidden for quintessence models with 
positive potentials, so a large fraction of early dark en- 
ergy is not allowed for quintessence in a flat universe. 

The same w = — 1 barrier for quintessence leads to pre- 
dictions of one-sided deviations from flat ACDM observ- 
ables when spatial curvature is allowed to vary. Closed 
universes reduce D*, so matching CMB constraints re- 
quires the ability to lower the dark energy density with 
w < — 1. Therefore, quintessence predictions from SN 
and CMB data consistent with w = — 1 favor open uni- 
verses. Curvature in open universes causes additional 
growth suppression, particularly at z < 2. 

Open quintessence models can also have significant 
early dark energy that suppresses growth by a constant 
factor at low redshift. Even with the additional freedom 
in curvature and early dark energy, SN and CMB data 
still provide general predictions for quintessence. For ex- 
ample, neither the growth relative to recombination nor 
absolute distances at low redshift can be significantly 
larger than in flat ACDM. Growth relative to the present 
out to z ~ 2 cannot differ from its standard behavior by 
more than ~ 2 — 4%. 

Allowing the low-rcdshift dark energy equation of state 



to vary beyond the range of quintessence enables new 
types of models since the absence of the w > — 1 bound 
permits deviations in dark energy density both above and 
below the constant density of flat ACDM. These general 
smooth dark energy models are therefore able to have sig- 
nificant early dark energy in a flat universe, and closed 
universes consistent with the SN and CMB data are also 
possible. As a result, smooth dark energy models can 
have both larger and smaller growth and absolute dis- 
tances relative to flat ACDM, unlike the one-sided pre- 
dictions of quintessence. 

Even the most general model class including large dark 
energy variations at low redshift, early dark energy, and 
nonzero curvature makes some generic predictions given 
future SN and CMB data. Growth at z < 2 normalized 
to its present value can be no more than ~ 5% larger and 
the expansion rate at z > 2 no more than ~ 5% smaller 
than their values in a flat ACDM cosmology. The growth 
index 7 as it is typically defined is tightly constrained 
by SN and CMB data in ACDM, but not for more gen- 
eral dark energy models since both time variation of the 
dark energy equation of state and the possibility of cross- 
ing Sl m (z) = 1 in closed universes weaken the relation 
between growth rate and matter fraction. Fortunately, 
growth measurements at different rcdshifts still give tight 
predictions since the growth evolution at low redshift de- 
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pcnds mainly on curvature. Additional tests outside the 
scope of this work, such as searching for scale depen- 
dence of linear growth, could also falsify all smooth dark 
energy models. Falsification of the most general smooth 
dark energy predictions would require new paradigms for 
cosmic acceleration and possibly even gravity itself. 
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APPENDIX A: FISHER MATRICES AND 
LIKELIHOOD FUNCTIONS 

In this Appendix, we describe the (future) cosmological 
data and priors that we assume for this study. We give 
expressions for the Fisher matrices for these data, which 
we use to compute the principal components of the dark 
energy equation of state, and for the likelihood functions 
that we use in MCMC analysis. 

The Fisher matrix for supernovae is [32| 



F, 



\ - _ 2 dm(z a ) dm(z a ) 
~ 2-^ Ua dO, d9, ' 



(Al) 



where m(z a ) = 5 log[Hodi,(z a )] + M is the average mag- 
nitude of the SNe in the redshift bin denoted by z a , 
a a is the error in the average magnitude, and M. 
\ 1 — 5 log(i?o/Mpc _1 ) + 25 is a constant related to the 
unknown absolute magnitude of the SNe. 

For the fiducial supernova data, we take the expected 
redshift distribution for SNAP [87| plus a low-z sample of 
300 SNe at 0.03 < z < 0.1. The SNAP magnitude errors 
include both statistical and systematic components: 



Az, 



0.15 2 



0.02- 



l + z 
2.7 



(A2) 



where N a is the number of SNe in each bin of width 
Az (Az = 0.1 except for the statistical uncertainties in 
the low-z SN bin, for which Az = 0.1 - zf^ in = 0.07), 
and Az su b is the width of the sub-bins used to smooth 
the distribution of SNe in redshift. We use 500 sub-bins 
up to z m ax = 1-7. The second term on the right hand 
side of Eq. (IA2[) models a systematic floor that increases 
linearly with z up to a maximum at z max of 0.02 mag per 
Az = 0.1 bin (l25j . 

For the Planck CMB constraint, we start with the 2x2 
covariance matrix C CMB for the parameters 



= {ln(ZVMpc),ft m /i 2 } 



(A3) 



Here is the comoving angular diameter distance to 
recombination. We ignore additional CMB information 



about dark energy from the ISW effect except in the cur- 
rent prior on early dark energy described below. The 
elements of the covariance matrix are C^ AB = (0.0018) 2 
and C*£, MB = (0.0011) 2 and Cf 2 MB = -(0.0014) 2 . Rotat- 
ing to the space of MCMC parameters, e.g. £>f u n [Eq. ijHJ)] 
gives F CMB = D[C' CMB ]- 1 L> T , where Dy = dfij/dOi. As 
we shall see, for the likelihood evaluation it is more conve- 
nient to project the MCMC parameters onto the original 
basis of Eq. (|A3[) . For a similar treat ment of CMB priors 
on dark energy models, see Refs. [82l Il09| |. 

Priors on additional parameters can be included by 
adding the assumed inverse covariance to the appropriate 
entry of the Fisher matrix. The Fisher matrix for the full 
set of parameters is 
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F, 



SN 



? CMB 

ij 
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(A4) 



The priors and the parameters in the Fisher matrix de- 
pend on the particular application and differ between the 
PC construction and the the likelihood analysis. The pro- 
cedure for computing PCs from the Fisher matrices and 
the assumed priors are described in Appendix [Bl 

For the MCMC analysis described in III CI we assume 
a Gaussian likelihood, C oc exp(— x 2 /2), described by 

(A5) 



X 



XSN 



XCMB 



which includes contributions from the SNAP SN data, 
Planck CMB data, and our external priors. 
We model the SN % 2 term as 

B 2 

Xsn= A -~c> ( A6 ) 
A = 5 y [Alog(gori L (za))] 2 

a < 

B _ 5 y A\og(H d L (z a )) ^ 

a a 

C = £^. 

a a 

where A\og(Hod^) refers to the difference between a 
model H d\j derived from MCMC parameters and the 
fiducial value. The variance er 2 is modeled as in Eq. (|A2|) . 



The B /C term in Xsn comcs from the marginaliza- 
tion over A4. Because of this marginalization, the SN 
data are insensitive to redshift-independent shifts in the 
magnitudes m(z a ) caused by changes in combinations 
of other cosmological parameters, e.g. f2 m and the PC 
amplitudes for w(z). This shift in magnitudes corre- 
sponds to multiplying the distances by a constant factor: 
H d h (z) (l + 5)H a d h (z). 

The CMB contribution to \ 2 is 



XCMB 



(A7) 



where 9i is the same as in Eq. (|A3|) and 66i = 9{ — 9i |fid is 
the difference in ml?* and fl m h 2 from the fiducial model 
with £>* derived from the MCMC parameters. 
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When computing \ 2 we se t the SN magnitudes, 
and fl m h 2 to have the exact values predicted for the fidu- 
cial cosmology without any scatter in the simulated mea- 
surements. The resulting constraints from the data are 
therefore expected to be centered on the fiducial param- 
eter values rather than shifted away from them by ~ 1 a. 
We do this because we are mainly interested in the width 
of predictions for observables and not their central values 
for a particular realization of the measurement errors. 

In the MCMC likelihood analysis we employ three ex- 
ternal priors based on current data to limit our study 
to reasonable cosmologies, and one internal theoretical 
prior: 
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where the terms on the right hand side respectively re- 
fer to an HST Key Pr oject prior on the Hubble constant 
of width a(h) = 0.08 102J, a BAO prior on the angular 
diameter distance to z = 0.35 with a{h\D(z = 0.35)) = 

0. 037, roughly corres pon ding to the constraint from the 
SDSS LRG sample |0|, and a WMAP prior on the 
fraction of early dark energy with <j(Qoe(z*)) = 0.025, 
base d on the constraints on early dark energy models in 
Ref. [lOll]. 

In addition to constraints from current data, the last 
term in Eq. (|A8[) includes theoretical limits on the dark 
energy equation of state, w m - m < w < w max - These lim- 
its are typically implemented as an infinite barrier in xt, 
corresponding to some range allowed by the model class, 

1. e. a top hat prior on the MCMC parameters. To com- 
pute these priors, we start with the projection of w(z) 
onto PC amplitudes, 



1 



N; 



z.PC 
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w(zj) - w&d]ei(zj), (A9) 



where Wfid = — 1 unless otherwise specified. By finding 
the values of w(zj) within the allowed range [u> m i n , «Wx] 
that maximize or minimize at, we obtain limits on the 
amplitude of each PC, a\~^ < on < aq , where 



a 



(±) 



1 



2N Z . PC 

^(^max 



E 



[(W min + ^■'max 
W m in)\ei(Zj)\]. 



2wfid)e l (zj) 
(A10) 



The width of this prior, a£ — a£ , depends on the width 
of the allowed range of w but not the value of Wfid- 

We find a second prior on {c^} using the fact that 
restricting the range of w(z) places an upper limit on 



wad 



From Eqs. © and ©, 
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The bounds on w(z) impose an upper limit on this sum, 
[w(zi) - w Rd } 2 < V max , (A12) 

Vmax = max[(w max - Wfid) 2 , (W min - Wfid) 2 ]- 
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FIG. 22: Examples of equations of state constructed from the 
15 PCs in Fig. (fl]), where the PC amplitudes satisfy the priors 
given by Eqs. (|A10|) and (|A13f) for -1 < w < 1 (solid black) 
and —5 < w < 3 (dashed red). 



Combining Eqs. (|A11|) and (|A12|) we find that the PC 
amplitudes must lie within a sphere in the parameter 
space: 



(A13) 



i=l 



i = l 



where N c < N z ^pq is the number of components in 
the truncated set of PCs used for likelihood analysis 
(see § III A|) . For a given allowed range of w, the con- 
straint on PC amplitudes from this inequality is strongest 
when the fiducial model lies in the center of that range, 
Wfid = (i« m m + w max )/2, but in general the limits in 
Eq. (|A10[) are the stronger of the two PC priors. 

Figure [2"2l shows examples of w(z) parametrized by 15 
PCs that satisfy the bounds in Eqs. (|A10|) and (|A13|) cor- 
responding to the quintessence model class, with w m i n = 
— 1 and u> max = 1) and the smooth dark energy class, 
with Wmin = —5 and u> ma x = 3. In each case, the sum 
of the 15 components is allowed to violate the bounds 
on if at some rcdshifts. For some models, the addition 
of higher-variance PCs can correct for these deviations 
from the allowed range so that w m i n < w < w max every- 
where. However, this is not necessarily true for all models 
that satisfy the PC priors; the most important property 
of these priors is that they retain any model that does 
obey the limits on w while excluding a large number of 
unacceptable models. For more details on the derivation 
of these priors, see Ref. [88j where similar bounds are 
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used in the context of reionization models to constrain 
the ionized fraction of hydrogen to the range < x e < 1 . 

When switching to alternate priors to test volume ef- 
fects related to curvature as described in $ 1111 CI we trans- 
form the PC priors either by adding an additional term 
to the likelihood in the MCMC analysis or by modifying 
the weights of samples in the chain as a post-processing 
step. In either case we convert the usual top-hat PC pri- 
ors to ones that are flat in the density of each PC at z max 
relative to z — [Eq. [TJ] by multiplying C or the sample 
weight by exp(— Xpc/2), where 



Xpc 



3a, 



dz 



1 



(A14) 



In addition to priors on the low-z equation of state, we 
place a top-hat prior on the early dark energy at z > z max 
corresponding to w m i n < Woo < min(0, io ma x)- Since the 
early dark energy MCMC parameter is actually e w °° , in 
practice the prior we use for MCMC is e -1 < e w °° < 
for quintessence and e -5 < e w °° < for smooth dark 
energy. 



APPENDIX B: PRINCIPAL COMPONENT 
METHODOLOGY 

To compute the principal components of w(z) from 
the Fisher matrices of Appendix [X] (with Ff™ 1 speci- 
fied below), we first take the A^ Zj pc X N z ^c submatrix 
of (i 7 ^*) -1 corresponding to the redshift binned dark 
energy equation of state, w(zi). We then invert the 
N z ,pc x Af Zi pc matrix to get F^j , which is the origi- 
nal Fisher matrix marginalized over everything except 
{w(zi)}. Finally, we compute the eigenvectors of FW, 
which are the PC functions and normalized as in Eq. (J3J, 
and the eigenvalues, which are the inverse variances of 
the PC amplitudes. 

We evaluate F^ at the same fiducial model as for the 
MCMC likelihood, i.e. flat ACDM with Q m = 0.24 and 
h = 0.73. This model is consistent with current data, 
which should minimize the number of principal compo- 
nents needed to accurately parametrize viable dark en- 
ergy models. The exact choice of fiducial cosmology is 
unimportant as principal component shapes do not vary 
greatly with changes in the fiducial model that are con- 
sistent with current data. 

Our default binning scheme is A^.pc = 500 bins be- 
tween z = and z max . These bins are fine enough to 
obtain reasonably continuous PC shapes and to allow 
varying w(z) at z < z^ n , which has important conse- 
quences as we describe later. 

We choose z max = 1-7 to match the assumed maximum 
redshift of the SN sample. A smaller choice of z max would 
not significantly change the PC shapes at lower z but 
would result in a less complete set of PCs due to neglect- 
ing some of the SN data. Increasing ,z max would require 
additional support from SNc or other data at z > 1.7 for 



the PCs to have any weight at higher redshift; the CMB 
distance constraint helps somewhat but is still only a 
single data point for constraining the additional w{zi). 
Furthermore, if the expansion history at low z is near 
the fiducial flat ACDM model then the lack of weight in 
the PCs at high z is mainly a consequence of dark energy 
becoming less significant as redshift increases (80j . 

Unlike the MCMC likelihood analysis, we do not in- 
clude external priors from current data or priors on PC 
amplitudes in the Ff^ lor term. However, we do use pri- 
ors that correspond to fixing certain parameters besides 
{w(zi}}, i.e. Q, m , f^K, and Woo. In the rest of this section 
we explain our choices of which of the other parameters 
are fixed and which are marginalized over. Note that 
we do not consider fl m h 2 here because it is nearly fixed 
automatically due to the constraint in F^ MB . 

The question of fixing or marginalizing parameters is 
essentially a question of whether or not we wish to in- 
clude degeneracies between those parameters and w(z) in 
the low-variance PCs that we retain for MCMC analysis. 
Marginalizing these other parameters eliminates from the 
low-variance PCs the modes in w(z) that have a degen- 
erate effect in the SN distances and CMB data but are 
not necessarily degenerate in other acceleration observ- 
ables. Since this marginahzation generally results in an 
incomplete PC basis, the better choice is typically to fix 
the non-w(z) parameters with F? Iim so that the modes of 
w(z) degenerate with them are assumed to be well mea- 
sured and therefore are included in the low-variance PCs. 
Fixing these additional parameters ensures that the PCs 
we use are as complete as possible in the acceleration 
observables, and so we fix fix and Woo when computing 
the PCs. However, n m is an exception to this rule where 
having a fully complete basis for w(z) is not desirable. 

There are two types of degeneracy that can in principle 
exist between f2 m and w(z) in the assumed SN and CMB 
data. First, the dark energy can mimic some fraction of 
the matter density (leaving Hod^ (z) u nchanged) by ap- 
proaching w = at high z (e.g. jl26l |). In the context 
of a spatially flat geometry with no early dark energy, 
this degeneracy is eliminated by the CMB constraints on 
and fl m h 2 . The impact of these constraints weakens 
if we allow for the freedom to adjust cither curvature or 
early dark energy. By constructing PCs with curvature 
and early dark energy fixed, we assume that the matter- 
mimicking mode of w(z) can be measured by the SN dis- 
tances and CMB constraints. This assumption ensures 
that the PCs are complete with respect to this mode, re- 
gardless of whether we marginalize or fix f2 m in the PC 
construction. This type of degeneracy between f2 m and 
w(z) is therefore properly included in the MCMC pre- 
dictions for model classes with curvature and early dark 
energy. 

A second type of degeneracy is introduced by the min- 
imum redshift for which SNe can be measured in the 
Hubble flow. The ability to determine H^d^z) from 
SN observations depends on how well we can anchor the 
relative distances to z = 0. As noted in Appendix [5] 
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marginalization over the nuisance parameter M. causes 
the SN observations to be insensitive to constant shifts 
in relative distances of the form 

H d L (z) -» (1 + (S)H 4W, (Bl) 
M -> M - 51og(l + 5). 

We are not free to change H^d^z) at all redshifts since 
Hodh{z) — z at low 2 independent of the cosmology. 
However, it is possible (if unlikely) that a large variation 
in w(z) near z = changes the SN distances by a nearly 
constant factor at all redshifts except at z < z^„, where 
is the minimum SN redshift. When including the 
CMB data this effect creates a degeneracy between w(z) 
and f2 m since the shift in Eq. (|Blj) requires changes in 
f2 m and i?o to satisfy CMB constraints on £>* and il m h 2 . 

If we fix f2 m when computing the PCs, then we assume 
that the behavior of w(z) at z < z^ n is well constrained 
and therefore large variations in the equation of state at 
low redshift are included in the low-variance PCs. In- 
stead, we choose to reduce the impact of this degenerate 
mode of w(z) by marginalizing Q m in the PC construc- 
tion. Our basis for w(z) is therefore incomplete with 
respect to this mode, but there are several reasons for 
neglecting large variations in w(z) at z < z^ n . One 
benefit of this approach is that by reducing the degen- 
eracy between f2 m and the PC amplitudes we improve 
convergence of MCMC samples; in the presence of the 
full degeneracy, it is difficult to obtain well-converged 
chains even for the simplest class of models with PC- 
parametrized w(z). Furthermore, apart from the Hub- 
ble constant itself and the interpretation of SN data as 
measuring H d^(z) as opposed to d^(z)/dL{z^ n ), accel- 
eration observables are not significantly affected by this 
degeneracy. Finally, despite marginalizing f2 m when com- 
puting PCs we still retain enough of this degeneracy that 
predictions for the z — ► behavior of D(z) and H(z) is 
appropriately uncertain and limited by our prior on Hq, 
as shown in Fig. [5l for example. 

Improving measurements of Hq beyond the current 
level would further limit the possibility of these ultra- 
low redshift changes in w. Conversely, in the absence of 
such variation in w(z), precision Hq measurements play 
the same role as low redshift D and H measurements in 
all cases considered in the main paper. 

In summary, when constructing the PCs we take priors 
that fix fix = and = — I but marginalize f2 m . We 
do not employ the additional current priors from BAO 
distance to z = 0.35, HST Key Project measurement of 
Hq, or WMAP limits on early dark energy that are added 
to the likelihood analysis. 



APPENDIX C: COMPLETENESS TESTS 

When making predictions for general classes of mod- 
els, we need to make sure that the parametrization we use 
has sufficient freedom to explore all types of effects that 



models can have on the acceleration observables. In this 
appendix we present several tests of the completeness of 
our parametrization. We begin by justifying the num- 
ber of principal components of w(z) used in the MCMC 
likelihood analysis. We then examine the sensitivity of 
our results to the choices of fiducial model and maximum 
redshift for principal components. Finally, we discuss the 
limitations of our early dark energy parametrization. 

The PCs form a complete basis for w(z) ordered by 
how well they can be measured by the fiducial SN and 
CMB data. This ordering allows us to truncate the set 
of PCs to some small number that have the greatest im- 
pact on the fiducial data. Retaining a limited number of 
PCs is a practical necessity to make parameter estimation 
feasible, but we must make sure that the higher-variance 
PCs that we ignore do not make significant contributions 
to the expansion or growth observables predicted in the 
main paper. 

It is important to emphasize that we do not expect 
or demand completeness in unobservable quantities like 
w(z) itself; the high-variance PCs that we neglect can 
have large effects on the equation of state, but since all 
of the observables contain integrals over w(z) the effects 
of these rapidly oscillating PCs (sec Fig. |T|) tend to cancel 
out for the redshift-dependent quantities of interest. This 
is especially true for distances and integrated growth, 
each of which involves essentially two integrals of w(z) 
over redshift. Completeness is more difficult to attain 
for observables with a single redshift integral such as the 
expansion rate and growth rate, but the practical require- 
ment of a large volume for such observations makes any 
rapid evolution with redshift unobservable in practice. 

Our basic strategy for determining the number of PCs 
required for completeness, N c , is to repeat the MCMC 
analysis for each class of models using varying numbers 
of PCs. As we increase the number of components of 
w(z), we expect the resulting predictions for observables 
to eventually converge once we have reached the neces- 
sary number of PCs. This approach makes the value of 
N c to some extent dependent on what we assume about 
the data and the allowed models. For example, it may 
be that the N c + 1 component can have a significant ef- 
fect on certain observables if its amplitude aN c+ i is un- 
constrained, but limits on this amplitude from the data 
and/or priors keep the impact of this eigenmode on the 
observables small. Similarly, it is possible that the value 
of N c changes as we add the additional freedom of spatial 
curvature and early dark energy to the baseline model. 

Since our definition of completeness is based on the 
precision of predictions for a variety of acceleration ob- 
servables, N c will generally differ from (and be larger 
than) the number of dark energy parameters that can be 
measured to some specified accuracy [13, HH, If 27} or the 
number of parameters required by the data in a Bayesian 
model selection sense [e.g.. [2ll [6l| . 

Figure[22]shows a comparison of observable predictions 
for the baseline quintessence class (Qk = 0, Woo = — f), 
using either 10 or 15 PCs. There is little difference 
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FIG. 23: Test of PC completeness for quintessence mod- 
els. Predictions for growth and expansion observables from 
MCMC with 10 (light gray) and 15 (dark blue) PCs for flat 
quintessence models ( — 1 < w < 1) with iCoo = — 1 (no early 
dark energy). 



between the two sets of predictions, suggesting that 
N c ~ 10 is sufficient for this model class. In contrast, 
predictions in this class using only 5 PCs are significantly 
tighter than those with 10 PCs. 

Including both curvature and early dark energy in the 
quintessence model class does not alter the agreement 
between the predictions for 10 and 15 PCs, as shown in 
Fig. [24] There is slightly more variation in the limits on 
observables at z > z max , but we do not expect perfect 
completeness at high redshift anyway due to our simplis- 
tic early dark energy parametrization (see below). Even 
in our most general model class where we weaken the 
quintessence prior to — 5 < w < 3 while continuing to 
include curvature and early dark energy, as in Fig. [231 
the predictions remain robust to increasing the number 
of PCs from 10 to 15. 

These comparisons indicate that N c ~ 10 is sufficient 
for completeness in all model classes that we study here. 
In the main sections of this paper we present results from 
the larger, "overcomplete" set of 15 PCs. 

In our predictions throughout this paper, we have as- 
sumed a particular flat ACDM model both for PC con- 



FIG. 24: Same as Fig. [23] for non-flat quintessence models 
with early dark energy. 



struction and for creating the fiducial SN and CMB data 
for MCMC likelihoods. Since the true cosmology could 
be somewhat different, we can ask how the predictions 
for observables would change had we assumed a different 
fiducial model. To test this dependence, we have redone 
the MCMC analysis using fiducial cosmologies with var- 
ious values of constant wad =/= — 1. For simulated SN 
and CMB data based on wra = —0.93 (approximately 
the 68% upper limit of combined constraints on constant 
w from WMAP and current BAO and SN data [HI), 
the predictions for growth and expansion observables in 
the context of quintessence models with either curvature 
or early dark energy have similar uncertainties to those 
with a fiducial ACDM cosmology (Figs. [6] and 0). The 
main effect of increasing w^d is that it slightly weakens 
the impact of the quintessence w — — 1 barrier; moving 
Wfid away from this barrier makes it possible to slightly 
reduce the dark energy density at z < z max from its fidu- 
cial level, and therefore the predictions for quintessence 
models include some features that were previously only 
allowed in the more general smooth dark energy class. 
For example, compared with the flat ACDM growth his- 
tory, flat quintessence models with early dark energy are 

allowed to have AG/G 5% at 68% CL instead of 

—2%, and curved quintessence models with no early dark 
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FIG. 25: Same as Fig. [23] for smooth dark energy models 
(— 5 < w < 3) including curvature and early dark energy. 



energy can have AG/G ~ 3% at z = (68% CL) where 
previously only downward variations in G(z = 0) were 
possible in this class of models. The ability to falsify cer- 
tain model classes therefore depends on how consistent 
future SN and CMB data sets are with the standard flat 
ACDM cosmology; any significant variation would be in- 
teresting in its own right and would make some changes 
to the model testing results presented here but not the 
methodology or the logic of the results themselves. 

Another technical issue related to completeness is 
whether our choice of z max = 1.7 affects the predic- 
tions for obscrvablcs. This choice enters into both the 
definition of w(z) principal components as discussed in 
Appendix [B] and the likelihood for MCMC as the maxi- 
mum redshift of the fiducial SN sample. To distinguish 
between the two, let us call the maximum redshift for 
PCs z^ x and for the MCMC likelihood z£ ax . Note that 
the choice of z max a l so influences our definition of "early 
dark energy" by setting the minimum redshift at which 



If we keep z^ ax = 1.7 but extend the PCs to z max = 2.5 
(assuming a flat SN distribution at 1.7 < z < 2.5 with 
the number per bin equal to the number at z = 1.7 in 
the original distribution), the resulting predictions for 
observables in the flat, no early dark energy quintessence 
class are similar to those in Fig. [5l The predictions are 
slightly weaker, particularly at 95% CL, due to the ex- 
tra freedom in w(z) at 1.7 < z < 2.5. The fact that 
width and redshift dependence of constraints on observ- 
ables change little with increased z max indicates that 
our results are not strongly influenced by the choice of 
^max = 1-7. For example, the tightening of constraints on 
G, D, and H at high z in Fig. [5] is more a consequence of 
the transition from accelerated expansion to deceleration 
at z ~ 1 than it is of setting z max = 1-7. 

As another test of sensitivity to z max , we set = 2.5 
as before and also extend the SN distribution for the 
MCMC analysis to z max = 2.5 (with constant number per 
bin at z > 1.7 as for the PCs). The resulting predictions 
for flat quintessence without early dark energy are nearly 
identical at z > 2 to the ones in Fig. [5] At lower red- 
shifts, observables are slightly better constrained due to 
the additional SN data; the largest changes are at z = 1, 
with new 68% limits of AG/G > -1% and AD/D < 1% 
(68% CL). Predictions for growth and expansion observ 
ables are relatively insensitive to increasing both z ) 
and z£„„. 



PC 
max 



At redshifts above our default choice of z max = 1.7, 
we make predictions for observables by specifying w(z > 
.Zmax) through an early dark energy ansatz with a con- 
stant equation of state, w = Woo. While Woo is not a com- 
plete parameterization of early dark energy, it does pro- 
vide guidance for predictions. For example, even though 
a constant equation of state at z > 1.7 is a poor fit 
to the Albrccht-Skordis model in which a quintessence 
scalar field has an ex pone ntial potential modified by a 
quadratic polynomial [128| . Woo as an effective parame- 
ter can nevertheless simultaneously fit the CMB distance 
to the required Planck precision and the growth function 
at z max to ~ 2% accuracy. The early dark energy pa- 
rameterization acts as a diagnostic for whether the dark 
energy can ever become a substantial fraction of the en- 



ergy density at z 



> 



given CMB constraints on the 



distance and energy densities at z ~ z*. Our philosophy 
is to use an incomplete but representative parameteriza- 
tion that can be used to monitor the need for early dark 
energy. If a substantial fraction of early dark energy is re- 
quired by observations under this parameterization, then 
more complete descriptions and more detailed observa- 
tions will be required. 
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